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Abstract. Preliminary results for the Dalitz plot distribution of 77' — > rjTtK decays in the frameworks of Large-/V f Chiral 
Perturbation Theory and Resonance Chiral Theory are given. We hope our results to be of some relevance for the present and 
forthcoming analysis of these decays at GAMS, CLEO, VES, KLOE-2, Crystal Ball, Crystal Barrel, WAS A, and BES-III. 
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INTRODUCTION 

The decays rj' — » r\nn are interesting for several reasons. First, due to the quantum numbers of the pseudoscalar 
mesons involved, if the decay proceeds through resonances these can be mostly of scalar nature. G-parity prevents 
vectors to contribute. Therefore, this decay is specially intended for the analysis of the properties of the /o(600) (or 
ct) resonance even though the ao(980) is also present and indeed dominant. Second, the presence of 77 and 77' in this 
reaction is ideal for studying the mixing properties of both mesons. Third, and more general, this decay allows to test 
Chiral Perturbation Theory (ChPT) and its possible extensions such as Large-iV c ChPT and Resonance Chiral Theory 
(RChT). For all that, precision measurements on 77 and 77' would be very helpful and provide useful information 
in our understanding of low energy QCD. In particular, their decays are perfectly suited to study symmetries and 
symmetry breakings in QCD and the simultaneous treatment of both the 17 and 17' imposes tighter constraints on 
theoretical approaches than just the 17. Accordingly, there is at present an intense activity in studying these processes. 
Recently, the GAMS-47T and VES Collaborations have measured the related Dalitz plot parameters (GAMS-47T for 
the 77' — > 777r°7r° channel [1] and VES for the 77' — > r\n + n~ one [2]) complementing older results reported by an 
early GAMS Coll. [3] and CLEO [4]. In the isospin limit the values of the Dalitz plot parameters should be the same, 
however the experimental measurements show some discrepancies among them. Then, there is a need of improvement 
and new measurements of theses parameters are foreseen at KLOE-2, Crystal Ball, Crystal Barrel and maybe WASA 
with improved statistics and a good understanding of the systematical errors. At KLOE the 77' is produced via the 
process e + e~ — > 4> followed by <t> —> rj'y. For the 77' — > r\n + n~ decay channel around 21K events are now on tape. 
The background is seen to be very low and these data can easily be used to determine the Dalitz plot parameters. 
The increased luminosity at KLOE-2 (a factor 3) will open new possibilities for 77' studies. For instance, a Monte 
Carlo simulation of the 77' — > 77 tt + ?t process shows that the detector has a good sensitivity to a [5]. At Crystal 
Ball, the 77' is photo produced, yp — > rj'p, and can be identified, e.g., from its decay 77' — > rjK°K° — > 67. About 10K 
events were expected during 2009 for the 77' — > i~in n Q process. At WASA, large samples of 77' will be produced 
in the reaction pp — > pprj' as soon as the new detector is developed. Instead, the reaction pd — > 3 He77' is not very 
suitable for 77' decay studies. Finally, about 27 million and 12 million decay events could be detected at BES-III 
each year for 77' — > r\n + n~ and 77' — > 77^°^°, respectively. On the theory side, the 77' — > r\%% decays have been 
studied within an effective chiral Lagrangian approach in which the lowest lying scalar mesons are combined into 
a possible nonet [7] and, more recently, within the framework of £7(3) chiral effective field theory in combination 
with a relativistic coupled-channels approach [8]. Here, we present tentative results for the Dalitz plot parameters of 
the decays 77' — s- r\nn in the frameworks of Large-7V C Chiral Perturbation Theory and Resonance Chiral Theory. This 
contribution is based on a more complete analysis that will be submitted in the near future [9] (see also Ref. [10]). 



1 Talk given at HADRON 2009: XIII International Conference on Hadron Spectroscopy, Florida, USA, 1 1/29-12/04/2009 



TABLE 1. Dalitz slope parameters (experiment and theory) for rj' — » Tjlfilfl (second and third columns) and rj' — > 
rjK + TZ~ (fourth and fifth columns), respectively. The isospin-symmetric values of Eq. (7) are shown in the last column. 



Parameter 


GAMS-4^ [1] 


Theory [8] 


VES [2] 


Theory [8] 


This work 


a 


-0.066±0.016±0.003 


-0.127 ±0.009 


-0.127±0.016±0.008 


-0.116±0.011 


-0.121 


b 


-0.063 ±0.028 ±0.004 


-0.049 ±0.036 


-0.106 ±0.028 ±0.014 


-0.042 ±0.034 


-0.015 


c 


-0.107 ±0.096 ±0.003 




+0.015±0.011 ±0.014 






d 


+0.018 ±0.078 ±0.006 


±0.011 ±0.021 


-0.082±0.017±0.008 


±0.010±0.019 


-0.075 



DALITZ PLOT PARAMETERIZATION 



The Dalitz plot distribution for the charged decay channel is described by the following two variables: 

X=^(T n+ -T«-), y=^L_JEJL_ 1 , (1) 

where T n ^ denote the kinetic energies of mesons in the rj' rest frame and Q = T 7] + T K + ± T K - = m^i — m n — 2m n . 
The squared absolute values of the two decay amplitudes (charged and neutral) are expanded around the center of the 
corresponding Dalitz plot in order to obtain the Dalitz slope parameters 2 [2]: 

\A(X,Y)\ 2 = \,^\ 2 [\+aY + bY 2 + cX + dX 2 ] , (2) 

where a,b,c and d are real parameters and \^V\ 2 is a normalization factor. For the charged channel odd terms in X are 
forbidden due to charge conjugation symmetry, while for the neutral c = from symmetry of the wave function. The 
Dalitz plot parameters may not be necessarily be the same for charged and neutral decay channels. However, in the 
isospin limit they should be the same. A second parameterization is the linear one [12]: 

\A(X,Y)\ 2 c< \ l + aY\ 2 + cX + dX 2 , (3) 

where a is a complex parameter. Comparison with the general fit gives a = 2Re(a) and b = Re 2 (a) +Im 2 (a). Both 
parameterizations are equivalent if b > a 2 /A. 



EXPERIMENTAL DATA 



The latest available experimental information on the Dalitz slope parameters is summarized in Table 1 . The analysis 
by the GAMS-47T Collaboration is based on approximately 15000 events [1]. Note that b is negative here and the fit 
is thus not compatible with a linear fit. If a fit is done with the linear parameterization one gets a = —0.042 ±0.008, 
which is in agreement with an early measurement from GAMS based on 5400 events that gave a = —0.058 ± 0.013 
(assuming Im(a) = and c = 0) [3]. Both analyses are for the rj' — > Tjn^n channel. Regarding the rj' — > Tjn + n~ 
decay, the result from CLEO based on 6700 events yields a = —0.021 ±0.025 (assuming Im(a) = 0, c = and d = 0) 
using the same linear fit as for the neutral decay channel [4]. The VES analysis is based on roughly 13600 events 
obtained from the charge-exchange reaction n~p — > rj'n and 6500 events from the diffractive-like production reaction 
K~N — > rj'n~N [2]. A fit using the combined data sets from VES gave the parameter values shown in Table 1. Again, 
b is found to be negative and thus incompatible with a linear fit, while c is consistent with zero. The previous work 
supersedes a first study [13], where the rjK + n~ Dalitz plot has been investigated with a sample of approximately 7000 
events obtaining a = -0.072 ±0.012 ±0.006 — this is the real part of a while Im(a) = 0.0 ±0.1 ±0.0 — in the 
linear and a = —0.120 ± 0.027 ± 0.015 in the nonlinear parametrization. The C-violation parameter was compatible 
with zero, c = 0.021 ± 0.024. In average, a = -0.059 ± 0.01 1 [12]. Finally, a measurement of a not included in the 
average is reported in Ref. [14], where a = —0.08 ±0.03 (assuming Im(a) = and c = 0) with about 1400 events. For 
the general parameterization, we see from Table 1 that there is some tension in all the parameters. Notice, however, 



2 The parameterization in Eq. (2) had been proposed in Ref. [11] with an extra term eXY. The analysis of Ref. [2] included this term in their fits 
and found that parameter e is consistent with zero. 



that the theory model results of Ref. [8] for the charged and neutral decay channels are compatible among themselves. 
While there is agreement in a between the VES fitted value and the result of Ref. [8] there is not such when compared 
to the GAMS-47C value. The same happens to d but this time the GAMS-47C and the theory model results agree, 
although are in conflict with the VES reported value. The case of b is less severe. The VES value is compatible with 
the others only at 2a. Finally, there is a variation between the measured values of c. Nevertheless, the large statistical 
error of the GAMS-47C result makes this statement not conclusive. 



In this contribution we study the decays 77' — > r\%% in two different chiral theories: Large-A^ Chiral Perturbation 
Theory (Large-jV e ChPT) and Resonance Chiral Theory (RChT). Both theories emerge from Chiral Perturbation 
Theory with the following modifications. ChPT is based on the spontaneous symmetry breaking of the SU(3)l x 
SU(3)r chiral symmetry exhibited by QCD at low energies down to SU(3)v- As a result, an octet of very light 
pseudoscalar bosons (n,K,r]$) with small masses compared to a typical hadronic scale such as M p appear. The 
pseudoscalar singlet 770 is not included in ChPT as an explicit degree of freedom since the U ( 1 ) axial anomaly prevents 
this state from being considered a ninth pseudo-Goldstone boson. However, in the large-A^ t limit the U(1)a anomaly 
is absent and the original chiral symmetry of QCD in the massless limit is enlarged to U(3)l x U(3)r. Large-A^ t ChPT 
is organized as a systematic expansion in powers of momenta, quark masses and 1 /N c — p 2 <~ m q <~ 1 /N c = 0(8) — 
describing the low energy interactions of n, K, 77 and 77' [15]. Then, in principle, it is perfectly suited for the study of 
77 and 77' hadronic decays. In the previous two frameworks lowest-mass resonances are not yet incorporated although 
it seems rather natural to expect an important impact of these resonances on the physics of the lightest pseudoscalar 
bosons. A systematic analysis of the role of resonances in the ChPT Lagrangian, the so-called RChT, was performed 
in Ref. [16]. In the following two sections we give tentative results for the Dalitz slope parameters predicted by these 
two theories. 



where F = F K at this order and is the 77-77' mixing angle in the quark-flavor basis. Taking F K = 92.3 MeV and 
the latest value = 40.4° from KLOE [17] one gets B(r\' -» riK + n~) ~ 10~ 2 which is very small compared to the 
experimental value B(r\' — > r]n + n~) = (44.6 ± 1.4) x 10~ 2 [12]. The reason for such a small value is the appearance 
of m 2 in Eq. (4), i.e. a chirally suppressed contribution. Moreover, the Dalitz plot distribution is constant in clear 
disagreement with the measurement of unambiguous slope parameters. At next-to-leading order (nlo) the 77' — > r]n + n~ 
amplitude involves the low-energy constants L5 , Lg and the combination 3L^ + L3 together with the OZI-rule violating 
parameter A2 [15]. The contributions from L^,L% and A2 are all proportional to m 2 and therefore chirally suppressed. 
The contribution from 3L^ + L3 is seen to be dominant and indeed, at this order, contains all the dynamical information 
encoded in the Dalitz parameters. For the numerical analysis, we use values for the pseudoscalar decay constants and 
mixing angles from Ref. [18] and the low-energy constants from Refs. [19, 20]. The values for the mixing parameters, 
which are extracted from phenomenology, are in reasonable agreement with the predictions of Large-jV e ChPT [21]. 
After expanding the exact expression of |A n i | 2 in terms of X and Y one finds the following relevant fact: |A n i | 2 only 
depends on X 2 in the m n ^Q limit (keeping m^ ^i ^ 0). This implies that terms proportional to Y carry extra powers 
of m n , thus making these terms numerically comparable to those of order X 2 . As a consequence, to be consistent 
our Dalitz expansion needs to include terms of order X 2 Y and X 4 denoted by K21 and K40, respectively. We also find 
combinations of slope parameters, such as a/d and JQo/Kai. which are independent of the low-energy constants and 
given as functions of the pseudoscalar masses alone. Preliminary values for the Dalitz parameters obtained from this 
Large-A^ t expansion are 



CHIRAL THEORIES 



Large-A/ C Chiral Perturbation Theory 




(4) 



a[Y] = -0.282 , d[X 2 ] = -0.083 , b[Y 2 ] - 10~ 3 , jc 2 i [X 2 Y] = 0.012 , k 4Q [X 4 ] - 10 



3 



(5) 



As seen, the value for d is in accord with the VES result while the values for a and b are not. This latter may be 
interpreted as an indication of the important role played by the a meson in this process whose contribution is encoded 



in the Y terms. Since the a pole is close to the allowed kinematical region a truncated expansion in Y would not be able 
to reproduce their effects. On the contrary, the parameter d encodes the contribution of the oq which is kinematically 
suppressed and thus well represented by the expansion in X. For those reasons, a better description of the role played 
by scalar resonances in 77' — > T]%% decays is desirable and this is the subject of RChT. 



Resonance Chiral Theory 

In this chiral theory the 77' — > r\%% decays are driven by the exchange of the a and fo in the s-channel and the 
do in the t- and M-channel. The contribution of scalar mesons in the RChT Lagrangian with a minimal number of 
resonance fields is represented by two operators whose coupling constants are cj and c m [16]. The largest contribution 
to the 77' — > 7771:71: amplitude comes from the terms, which are proportional to the external momenta, while the rest 
is proportional to m 2 and then suppressed. At low energies, when s,t ,u <C Mf, the amplitude in the RChT framework 
recovers that of Large-A^ t ChPT up to subleading contributions in 1 /N c . After using the low-energy constants relations 
3L2 +L3 = c^/2M|, L5 = CrfCm/M^ and L% = c^,/2M|, and the large-iV c result A2 = 0, the two amplitudes match. For 
the numerical analysis, we use the same set of input values as for the Large-A^ ChPT case together with Ms = 980 
MeV for the scalar multiplet mass. We also use = 32 ± 13 MeV and c m = 42 ±7 MeV from Ref. [16]. They are 
obtained from the phenomenological values L r 5 (M p ) = (1.4±0.5) • 1CT 3 and L\(M p ) = (0.9 ±0.3) • 10~ 3 assuming 
scalar resonance saturation. These values of and c m allow for a prediction of B(rj' — > r\n + n~) in agreement with 
experiment. However, since the total decay width is largely proportional to c A d the propagated error is huge. Therefore, 
we prefer to fix cj from the experimental branching ratio [12], thus obtaining c L i = 28.9 ±0.2 MeV, and c m from 
the short-distance constraint in the single-resonance approximation Acd.c m = f 2 ~ f% [22]. The values for the Dalitz 
parameters one gets in this case are 

a[Y} = -0.116, d[X 2 } = -0.054 , b[Y 2 } ~ 10~ 3 , jc 2 i [X 2 Y] = -0.006 , k 4Q [X 4 ] ~ 10~ 3 . (6) 

Now, the values of a and d agree with the VES results while b is still very small. This behavior may indicate that a 
meson effects are already accounted for but not yet in precise form. So far, subleading effects in large-A^ such as final 
state interactions (making the a a broad resonance) or the use of different scalar masses within the multiplet are not 
considered. Different ways of including these effects haven been studied in Ref. [9]. For instance, the inclusion of a 
a-fo mass difference along with a self-energy in the a propagator that reproduces the correct a pole in the complex 
plane (the same is done for oq and /o) gives 

a[Y] = -0.121 , d[X 2 } = -0.075 , b[Y 2 } = -0.015 , jc 2 i [X 2 Y] = -0.008 , k 4Q [X 4 ] - 10~ 3 . (7) 

We see that a and d are in nice agreement with the VES results, whereas the prediction for b, although increased in 
one order of magnitude, is still far from its experimental value. For completeness, we display in Table 1 our former 
results for the parameters a, b and d. 

SUMMARY AND DISCUSSION 

In this contribution we have computed preliminary values for the Dalitz slope parameters related to 77' — > 777^ decays 
in the frameworks of Large-A^ c Chiral Perturbation Theory, at lowest and next-to-leading orders, and Resonance Chiral 
Theory in the leading 1 /N c approximation. We have seen that next-to-leading order Large-A^ t ChPT fails to reproduce 
the Dalitz parameters associated with scalar exchanges in the ^-channel, thus indicating the relevance of a meson 
effects in this channel. The contributions from the t- and M-channel are, however, well reproduced. The prediction for 
the Dalitz distribution is improved in the context of RChT where scalar resonance effects are included explicitly and 
not integrated out in the low-energy constants. We have then seen that the contribution of the ao is dominant but that of 
the a is essential to reach agreement at the current precision level. The situation gets even better as soon as subleading 
effects in 1 /N c are considered. The incorporation of a a-fo mass splitting and some final state interactions through 
the modification of the a propagator seem to be unavoidable in order to match the a and d parameters of the Dalitz 
plot expansion. The b parameter is, however, not yet reproduced. For the time being, we do not foresee how to make 
this parameter to agree. This is in contrast with the predictions of the 1/(3) chiral effective field theory with coupled 
channels where a and b are nicely reproduced but d is not. Likewise, current experimental results show some tension 



among themselves, especially in the case of the a parameter. We hope that forthcoming analysis of the 77' — > T]%% 
decays at different facilities will clarify the present situation and allow for a distinction between different chiral theory 
predictions. 
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